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Abstract. Starting from the exact non- linear description of matter and radiation, a 
fully covariant and gauge - invariant formula for the observed temperature anisotropy of 
the cosmic microwave background (CBR) radiation, expressed in terms of the electric 
(Eab) and magnetic {Hab) parts of the Weyl tensor, is obtained by integrating photon 
geodesies from last scattering to the point of observation today. This improves and 
extends earlier work by Russ et al where a similar formula was obtained by taking first 
order variations of the redshift. In the case of scalar (density) perturbations, Eab is 
related to the harmonic components of the gravitational potential $fc and the usual 
dominant Sachs -Wolfe contribution ^T^/Tr ~ <&fe to the temperature anisotropy is 
recovered, together with contributions due to the time variation of the potential (Rees - 
Sciama effect), entropy and velocity perturbations at last scattering and a pressure 
suppression term important in low density universes. We also explicitly demonstrate 
the validity of assuming that the perturbations are adiabatic at decoupling and show 
that if the surface of last scattering is correctly placed and the background universe 
model is taken to be a flat dust dominated Friedmann- Robertson- Walker model 
(FRW), then the large scale temperature anisotropy can be interpreted as being due to 
the motion of the matter relative to the surface of constant temperature which defines 
the surface of last scattering on those scales. 



1. Introduction 

The study of the cosmic microwave background (CMB) radiation is the corner stone of 
modern Big Bang cosmology and has led to its widespread acceptance. Over the last few 
years improved measurements of the temperature spectrum and anisotropy has led to a 
better understanding of the origin and evolution of large scale structure in the universe. 

The physical basis for using anisotropies in the CMB to constrain competing theories 
of galaxy formation is the Sachs -Wolfe effect pH]. The basic assumption is that the 
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observed CMB photons travel to us without significant interaction with matter from a 
redshift of about z = 1200. At redshifts greater than this, the universe is ionized and 
photons are coupled to the electron - baryon plasma through Thompson scattering. 

The process of decoupling i.e. the transition of the CMB from a coUisional regime to 
being free photons does not take place instantaneously. The thickness of the decoupling 
shell Az is approximately 1/15 of the mean redshift, which from our point of view as 
observers is relatively narrow so for many purposes we can treat this shell as a 
sharp surface, called the surface of last scattering (SLS). The correct way of placing 
this surface is by determining where the optical depth due to Thompson scattering is 
unity [|^, |10|- This occurs, to first order, where the radiation temperature 

(which is equal to the matter temperature in the strongly coupled region prior to 
decoupling) reaches the matter ionization temperature, so if we take decoupling as 
happening essentially instantaneously, the SLS is, to good approximation, a surface of 
constant radiation temperature. 

Causally connected regions at the SLS, as viewed by an observer today, subtend an 
angle 9 ~ 2^/^^^, where Qq is the present value of the density parameter, so large scale 
anisotropies on angular scales greater than 7 degrees are unaffected by the small scale 
physics of decoupling, and so represent primordial perturbations. These anisotropies 
arise because photons traveling from the SLS are red -shifted slightly more than they 
would be in a perfectly homogeneous universe as a result of having to climb from an 
increased gravitational potential due to density perturbations over the surface. 

The calculation of CMB anisotropies on angular scales larger than a degree, 
corresponding to scales larger than the Hubble radius at the time of decoupling, is very 
simple in principle, however in practice the formulation of a completely self - consistent 
theoretical picture has been plagued by a misunderstanding of the observational meaning 
of certain temperature perturbation measures and whether or not these measures are 
gauge -invariant. Furthermore many approximations are made in deriving the various 
contributions to the anisotropy, without adequate consideration of whether or not these 
approximations are justifiable or consistent with each other. 

In this paper we attempt to address some of the above problems, by deriving from 
first principles a formula for the CMB anisotropy which is both simpler and easier to 
interpret than many of the usual treatments. The approach we take is to integrate 
photon geodesies from the time of last scattering to today, obtaining a general (model 
independent) formula for the observed temperature in a given direction in the sky and 
expressing the result in terms of covariantly defined gauge - invariant quantities. The 
temperature anisotropy is then given by subtracting this result for two independent 
directions corresponding to different points of emission on the SLS. This improves and 
extends earlier work by Russ et al |]2^ by (i) starting from a general non- linear treatment 
of the geodesic equation and linearizing to obtain the almost FRW result, rather than 
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taking first order variations of the redshift; (ii) correctly defining the temperature 
perturbation measure 5T /T] and (iii) performing a two - fiuid analysis instead of treating 
the radiation as a test field on a single -fiuid background. 

Following this, we demonstrate that given a number of "standard" assumptions, the 
classical Sachs -Wolfe result is recovered and we also consider in detail the validity of 
the assumption that the perturbations in the total energy density of the photon - baryon 
fiuid at decoupling are adiabatic. 

Conventions on signature, Riemann and Ricci tensors are as in |]^, and the speed 
of light is taken to be unity (c = 1). Standard General Relativity is assumed, with 
Einstein's equations in the form Gab = i^Tab where Gat is the usual Einstein tensor, 
K, = 8ttG is the gravitational constant and Tab is the energy momentum tensor of the 
matter. Most of the notation is the same as in P, Q and any changes are stated in the 
text. 



2. Basic equations and notation 

For the sake of self - consistency and clarity, we will first summarize the covariant 
approach to cosmology. 



2.1. The covariant approach 

As in Hawking |]14| and Ellis [§, |^, the hydrodynamic fiuid 4- velocity (tangent to 
the worldlines of fundamental observers in the universe) is = dx°'/dt {u°'Ua = — 1), 
where t is the proper time along the fiow lines. The projection tensor into the tangent 
three - spaces orthogonal to m° (the local rest frame of these observers) is 

hab = 9ab + UaUb ■ (l) 

The first covariant derivative of Ua can be uniquely decomposed into four parts: 

Ua;b = ^^^VfoMa " UaUb , (2) 

where 

^^VfoMa = (Tab + i^ab + \Qhab , (3) 

and '■^-'Va is the spatial gradient operator, orthogonal to u'^. Here 9 = M°;a is the 
volume expansion, aab = cr^ab) is the shear tensor {aabu"" = cr^a = 0), Uab = ^^[ab] is 
the vorticity {ujabu^ = 0) and = Ua-bu'' is the acceleration (the dot denotes a proper 
time derivative). It is useful to introduce a length scale along the fiuid fiow lines by the 
relation 
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When the universe is an exact FRW spacetime, H is just the usual Hubble parameter. 
In general, however, the evolution equation for the expansion B is the Raychaudhuri 
equation 

e + + 2((t2 - Cj2) - (3) V^Ma - Uail" + + 3p) = , (5) 

where = \aabC^°''^ and u"^ = ^ojab^""^ are the shear and vorticity magnitudes and /z and 
p are the energy density and pressure respectively. 



2.2. Matter and radiation 

Fixing u"" so that it corresponds to the Landau- Lifshitz (energy) frame |l^t and 
considering only small deviations from equilibrium (so that the velocities of the matter 
and radiation relative to this frame are small), the total energy momentum tensor for 
matter (m) and radiation (r) is given by 

Tab = fJ'UaUb + phab + VT^fo'' , (6) 

where 

Ai = At(m) + /i(r) , P = |/i(r) , (7) 

(r) 

and Tc^f, is the anisotropic pressure of the radiation. 



2. 3. Component fluid equations 

Relative to this frame the conservation of energy and momentum for non- interacting 
matter and radiation are given by 

+ |/i(,)e + (3) V'^gi'-) + 2q(;)u- + 7riJa«^ = , (8) 



and 



AM + /iMe + (3)V'^gf) = 0, (9) 

h'^aqi'^ + l^i(r)Ua + l^'^Vafi^r) + V'vri^ (10) 
+ U\^f! + [aab + Wab + l&hab) = , (11) 



where 



h\qi"''^ + ^{m)Ua + (cTab + OJab + |B/lafe) q\m) = , (12) 



and and l/j'") are the velocities of the matter and radiation relative to m": 

= - Ua , Vi™) = - Ua . (14) 
f In this frame the total energy flux vanishes. 



2.4- Total fluid equations 

Because we have chosen to work in the energy frame, the conservation equations for the 
total fluid are considerably simpler than those for the individual components: 

fi + hQ + TrlJfx'^' = , (15) 

hua + ('V.p + VVlJ + ^^1? = , (16) 

where 

h = fJ'{m) + |/^(r) , (17) 

is the sum of the total energy density and pressure. 
2.5. FRW models 

In the case of a FRW universe, m^"^) = -u"^) = m° and vr^j, = 0, so the energy momentum 
tensor (|]) necessarily reduces to the perfect fluid form 

Tab = f^UaUb + phab, (18) 

and 

Ua = ^'^^aP = , (19) 

SO dynamics of matter - radiation models are completely determined by the energy 
conservation equations 

+ 4/i(,)e = , (20) 

A(m) + 3/i(m)6 = , (21) 
/i + /iB = , (22) 



H^ + — = lKfx, (23) 



together with the Friedmann equations 
K 

where K is the spatial curvature constant, and 

3H + 3H'^ + Ik{i2 + p) = , (24) 
which is (B) specialized to a FRW model. 
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Figure 1. A diagram illustrating the geometry of the Sachs- Wolfe effect. It clearly 
shows that the observed temperature anisotropy at xr is given by the difference in 
temperature of the CMB observed in different directions on the plane of the sky, 
corresponding to different points on the surface of last scattering (SLS). 



3. A gauge - invariant measure of CMB temperature anisotropies 

The aim of this paper is to derive in the simplest possible way a covariant and gauge - 
invariant formula that describes the observed CMB temperature anisotropy ATr/Tr. 
This improves on earlier work by Russ et al where a similar formula was obtained 
by taking first order variations of the redshift. 

Before sketching this derivation, let us first clarify the difference between the gauge - 
invariant temperature perturbation 5Tr/Tr and the observed temperature anisotropy, 
since this continues to be a source of considerable confusion in the literature. 

If the point of observation is defined by the spacetime point x^, with coordinates 
x^, then the temperature measured by an observer at x/j in a direction e^, with 
components e% can be decomposed into the average bolometric temperature Tr(x^) 
and the covariant and gauge - invariant temperature (5TR(x",e^) [p!9[] : 



Tr(xr, e^) = rR(x/j) + 6Tr{:s.r, e^) , (25) 
fnix'},) =^ [ Tnix%e'\)dn. (26) 

47r J An 



The observed temperature anisotropy is given by 
ATji Tr{xr, Ba) - Tn{xR, e^) 



(27) 



where e^ and e^ correspond to two different directions on the plane of the sky 
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corresponding to two different points of emission and on the last scattering 
surface (see figure 1 above). If 6Tr <C 1 this becomes 

ATr 5Tr 5Tr 

-7fr- = -j^i^R, ga) - -j^i^R, gb) , (28) 

-'■ R -'■ R R 

since the average bolometric temperature contributions to T cancel out when subtracted 
at the point of observation x^. 

4. A gauge - invariant formula for 5Tr/Tr 

Having clarified the meaning of temperature perturbations and the observed 
temperature anisotropy, we will now proceed to derive from first principles an expression 
for STr/Tr in terms of the covariant quantities defined in section 

In principle what one should do is integrate the Louville equation in curved 



spacetime for a generalized gauge - invariant distribution function ||28|, |T9[ first through 
the decoupling phase and then from decoupling until today to obtain the perturbed 
spectrum of photon energies observed in a given direction on the sky. The 
temperature anisotropy is then obtained by subtracting these results for two independent 
directions. 

It turns out however, that because photons are essentially coUisionless after last 
scattering, the complete physical content of the Louville equation is contained in the 
geodesic equation, so it is much simpler instead to integrate the photon energies E up 
null geodesies connecting points of emission on the last scattering surface with the point 
of observation x/j here and now. 

Photons moves on null geodesies x°'{X) where A is an affine parameter. The null 
vector to these geodesies is 

P^ = ^^ P^Pa = , (29) 
and satisfies the geodesic equation 

= . (30) 
The 3 + 1 decomposition of this null vector 

p'' = E{u^ + e'') , e'^u, = 0, e-^e, = 1 , (31) 
defines the energy 

E = -p^Ua (32) 



of a photon relative to the four - velocity of the observer and its relative (spatial) 
direction of motion e"". 
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The temperature Tr = TR(x^,e£;) of the CMB observed at the reception point R 
with spacetime coordinates in a given direction Ge is given by 

Ta 1 + Z {p''Ua)A ' ^ ^ 

where Ta = Ta(xa) is the temperature at the emission point and z = z{eA) is the 
redshift between emission and reception. 

The equation describing the variation of a photon's energy along a null geodesic 
(parameterized by A) is obtained by differentiating equation ( p2D with respect to A and 
using the geodesic equation (|30D. This yields 

" ' (34) 



dX 



(35) 



Substituting for Ua-^ from equations (0) and (^ gives 
dX 

At this point one could already integrate for E, however for the problem that we wish 
to discuss, it is more convenient to substitute for the expansion in terms of known 
gauge - invariant perturbation variables. This is best done by projecting the spatial 
gradient of the radiation energy density Xj^^^ = '•^•*Va/i(r) along the null vector p": 



(36) 



dx- ' dX 
Substituting for from the energy conservation equation for the radiation (|2.3| ) one 
obtains 



P 



dX 



It follows that 



(37) 



, (38) 



Mr) I 

and substituting this into equation ( P5| ) gives a completely general equation for the 
variation of a photons' energy along a null geodesic: 

1 dE 1 



E dX 4/i(r) dX 
where E is given by 



(39) 



4/i(r) 
1 



E 

4/i(r) 



+ -^(^abP^P^ + UaP" ■ 



(40) 
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Integrating up a null geodesic from the point of emission at last scattering to the 
point of reception and using equation (|33D one obtains an exact formula for the 
temperature at reception Tr: 



In 



Tr 
Ta 



In 



FdX . 



(41) 



Substituting for Tr from equation (25) and using the Stephan - Boltzmann law 
aT"^ we obtain 



In 



1 + 



R 



T, 



R 



R 



FdX . 



(42) 



It is important to realize that apart from assuming that photons are collision free, no 
approximations have been made in this section so far, and the above result is therefore 
valid for any choice of background geometry and matter descriptionf. 



Linearization about FRW models 

To obtain an expression for 5Tr/Tr valid in an almost FRW model (a spacetime where 
these variables are small) we approach this universe from the equations valid in a general 
spacetime rather that adopting the standard procedure of perturbing an exact FRW 
model. The linearization procedure we apply consists of dropping terms such as T^ab C 
in equation (40), i.e. terms which are second order in the gauge - invariant variables, 
retaining only those terms which are linear, for example X^^^ Linearizing equation 
(1^) following this procedure yields the following result 

STn ( 1 . 1 



where P^^^ = aX^^ j and we have used the normalized background (FRW) expression 
for the photon energy E = 1/a. 

This formula expresses the generation of CMB anisotropies by cosmological 
perturbations in it's clearest form with each term having a direct physical interpretation. 
First it should be noted that T)^'', V^^'^ and iia contain both a scalar and vector 
part, while the shear aab is made up of contributions due to scalar, vector and tensor 
perturbations. Focusing on scalar perturbations, V'^^ and V^"^^ characterize density and 
velocity perturbations in the radiation relative to u"', the acceleration term u"" represents 



possible pressure suppression effects ||2^ and the shear relates to perturbations in the 
gravitational potential. 

f In the case of a FRW model, F = and ^(r) a"^ so the standard result of Tr/Ta — aAjo-R is 
recovered. 
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We can express this result in terms of total matter variables by using the following 
results 0: 

+ = ^ (i - + l^s^r' , (44) 



y(r) ^ [H^^yirm) ^ ^45^ 



where 



cl = (46) 

3(4/i(r) + 3/i(„)) 

is the speed of sound in the total fluid, 

f,V, = + fi^r)Vj[^ (47) 

is the total perturbation in the energy density and 

gir^) = - , l/i"""^ = - u'r^ (48) 

d ^ a ^ a ' a a a \ / 

are the entropy and relative velocity perturbations respectively @]. 

At the time of decoupling, if the present value of the density parameter Qq > 0.1, 
the universe is matter dominated to a good approximation, so h — > l^{m) 

and cl 0. It 

follows that above results reduce to 

+ au, = h), + \s^r^ , Vi^^ = Vi^-) (49) 



and the expression for 5Tr/Tr becomes 

^ A- { —V + aaabpy' ] d\ , (50) 



J-R 

where 



A \3a 



A= - r ^ (^(^'"V" - (3) v'^vj™)) d\ . (51) 



I A 3a 

We thus have two contributions: one due to perturbations in the total energy density 
and pressure and the other. A, due to the difference in the dynamical behavior of the 
matter and radiation density and velocity perturbations. 



5. The temperature anisotropy due to gravitational potential 
perturbations 

In this section we will deal with the contribution to CMB anisotropics due to 
gravitational potential fluctuations. In order to do this it is first convenient to write the 
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formula (|50| ) in terms of the electric Eat and magnetic Hat parts of the Weyl tensor. 
This is done by using the two linearized Bianchi identities which relate to Eab 0: 



(52) 



(53) 



to substitute for the shear and Va in (|50D. This leads straightforwardly to the following 
result: 

5Tr . fR 1 



A 



V^Eabp' - 2a + 3HE,b) P^'p" 
dX . 



(54) 



This expression is closely related the formulae derived by Magueijo [20| and Durrer 



In the case of scalar gravitational potential fluctuations the magnetic part of the 
Weyl tensor Hab vanishes (since it only contributes to vector and tensor perturbations 
0, |13[), so (54) reduces to 

^-^ = A-r^ l^'^V'Eabp' - 2a (Eab + SHEab) dX . (55) 

For scalar perturbations, the electric part of the Weyl tensor (with wave number fc) is 
related to the harmonic component of the perturbed gravitational potential $fc = $fc(t) 
as follows HI: 



Eab = -^^kQab 



(56) 



where Qab is a covariantly defined scalar harmonic and k is the eigenvalue associated 
with Qab (it is the wave number if the background FRW spacetime is flat). A discussion 
of these harmonics and their properties is given in appendix A and . 

Substituting ( |5BD into (^) and using equations ( |A3| ) and ( [A5| ) in appendix A, we 
obtain 



--A-2 



A a-'[l{m) 



(a$fcQ)' 



K{a<!>k)Q' -a^ (a$fc)' (aQ 



dX 



(57) 



where the prime denotes differentiation with respect to A. In the background FRW 
model the energy conservation equation (PTj) can be integrated to give 



aa 



(5^ 



where fiE and are the background values for the energy density and scale factor at 
the time of emission. 
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Substituting (|^) into (57) and integrating the first term by parts gives 

+ - r \kQ + (aQ')1 {h^u + d\ , (59) 

where we have dropped the term evaluated at reception since it has no angular 
dependence. This result is true for a general FRW background. The $fc represents the 
integrated or Rees- Sciama effect which is important if the potential is non- stationary, 
for example in open universe models. 

5.1. Large scale temperature anisotropies 

To further simplify the above problem, we will now make the standard assumptions of 
assuming that the background is a flat {K = 0) FRW model, and consider temperature 
anisotropies arising as a result of density perturbations on scales much larger than 
the Hubble radius at decoupling. In this case the Friedmann and energy conservation 
equations (^- 0) lead to the following background evolution for the scale factor: 

a=(/3t)i , /?' = fa, (60) 

and in the matter dominated regime the potential fluctuations $fe satisfy the following 
differential equation [|l^: 

^k + ^H^k = 0, (61) 

which follows from ( ^2]) and the shear propagation equation when Hat = 0. Substituting 
for the scale factor (|60|) in (59) and integrating by parts, using (|6TD to substitute for 
second derivatives in gives the following result: 



+ 1^^ (f) ' {a^kQaP"")^ + lA^kQdt . (62) 



where 



is the ratio of the Hubble scale Xh = 1/H io the comoving scale A at the time of 
decoupling and the last term represents the integrated Sachs- Wolfe effect. 

If the gravitational potential is approximately constant and we consider scales 
much larger than the Hubble radius at decoupling, so that TIa -C 1, we recover the well 
known result of Sachs and Wolfe, together with additional contributions due to velocity 
and entropy perturbations 

^-^ = A+{^uQ)a- (64) 

J-R 
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Finally if we subtract this result for two independent directions {A and B), the observed 
temperature anisotropy ATji/Tji is obtained: 

AT 

-^ = AA+A{<l>kQ), (65) 

where 

A ($feg) = (<i>,g)^ - {<^kQ)B ■ (66) 

is the difference in the gravitational potential between separate points A and B on the 
SLS. 



6. The adiabatic assumption 

One of the most common assumptions made when discussing large scale CMB 
anisotropics is that the perturbations in the total energy density are adiabatic at the 
time of decoupling. Let us now consider whether or not this assumption is consistent 
and how it affects the results presented above. In order to achieve clarity on this issue, 
we need to consider (i) the large scale evolution of density, entropy and relative velocity 
perturbations during the collision dominated period prior to decoupling as these provide 
the initial conditions at the time of decoupling; (ii) how these initial conditions relate 
to the correct placing of the surface of last scattering and (iii) the evolution of these 
perturbations after decoupling. 



6.1. Before decoupling 



On large scales, in the matter dominated limit, the dynamics of density A(j.) = 
a(3)v«r'^'-), entropy S(^rm) = a'-^^VS^^'^^ and relative velocity V(rm) = a^^) V^l/j^™) 
perturbations in a photon - baryon universe are described by the following set of 
equations M: 



A(r) + 2HA 



(r) 



-hA 
2 



(r) 



(rm) 



H I 



h 



■Rr S, 



(rm) 



(67) 



and 



Sirm)+H\- — 



-Rc + 1 'S'l',., 



(rm) 



(rm) 



3 h 



+ 



-Rc V( 



1 1 



(rm) 



A 



(r) , 



(68) 



(69) 



where Rc{t) = 1/Htc is the ratio of the horizon size to the mean free path for photons 
colliding with electrons and Tc is the mean collision time of photons with electrons. 
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Equations (0-|69|) can be integrated to give solutions for density, entropy and 
velocity perturbations in the tightly coupled regime before decoupling: 

A(,) = t'/'Aa , (70) 

S(^rm) = So + B f"^ -e-^(*)rft , (71) 
Jo a 



V(^rm) = C7(t)e-^W . (72) 

Before decoupling Rc{t) is much greater than unity since the mean collision time between 
photons and baryons tends to zero, and therefore P{t) ^ 1 in this limit. This means 
that the solution for S(^rm) (|7l|) settles down to a constant value immediately after it 
is provoked, and so entropy perturbations have essentially one mode which is constant 
in time. This is due to the fact that the matter and radiation are so tightly coupled 
that the matter cannot move relative to the radiation. This behavior can be seen by 
looking at the solution for V(rm) (172]) which is exponentially driven to zero. These 
solutions imply that if the perturbations are initially adiabatic, as suggested by many 
inflationary scenarios, they will remain so until the time of decoupling. 

6.2. Adiabatic perturbations at decoupling 

Given that the pre - decoupling perturbation dynamics can lead to adiabatic initial 
conditions, let us consider whether they are compatible with the proper placing of the 



SLS |2g, |2|, |T0| 



The correct way of placing this surface is by determining where the optical depth 
due to Thompson scattering is unity. This occurs, to first order, where the radiation 
temperature, which is equal to the matter temperature in the strongly coupled region 
prior to decoupling: T(r) = T(m) = T, reaches the matter ionization temperature, so the 
last scattering event A on each null geodesic is characterized by 

Ta = Tion ■ (73) 

Thus if we take decoupling as happening essentially instantaneous, the SLS is, to good 
approximation, a surface of constant radiation temperature and so by the Stefan - 
Boltzmann law /i(^) = aT^, also one of constant radiation density: 

AT = ^ ^'i"^ = , (74) 

where AT = Ta—Tb is the difference in temperature between separate points of emission 
xa and xb on the SLS (see figure 1) and V'^J'^ is the spatial variation of {X(^r) orthogonal to 
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the normals to the surfaces of constant radiation density (V^'^n"' = 0). Transforming 
to the energy frame u^, we can relate V^J'^ to V^J'^: 

Vj[^ = + AaHVa , V'' = u''-n\ (75) 

and taking its spatial divergence we obtain the corresponding result for scalar 
perturbations: 

A(^) = A(^) + AaHV , (76) 

where 

A(,) = a(3)V"Pi") , V = a^^^V'Va . (77) 

If the perturbations are adiabatic at decoupling: 

4 

S{rm) = ^ A(r) = -A(m) , V(^rm) = , (78) 

SO equation (^) becomes 

A(,) = ^A(^) + AaHV . (79) 

Hence, for a surface of constant radiation density A(^.) = 0, which defines the SLS, we 
find that: 

A(^) = -3HaV . (80) 
Using equations ( ^31) and we can relate A(m) to the electric part of the Weyl tensor: 

^2{3)ya(3)y6^^^ = |K/iA(„) , (81) 

and combining this with ( pO]) we obtain: 

(3)ya(3)yb^^^ = ~ifV . (82) 

Using ( p6D and the results in Appendix A, the LHS of (^) can be written in terms of 
the harmonic components $fc of the perturbed gravitational potential, while the RHS 
can be decomposed into its harmonic components Vk (see equation 102 in 0]): 

(3)v-(3)v^E,, = 15 ($fcQ) , V = -kVkQ . (83) 

Substituting these results into ( [8^ ) and using equation ( [63| ) gives the following result: 



^kQ = \n-^VkQ . (84) 

It therefore follows that for adiabatic perturbations the large-scale temperature 
anisotropy 

^ = A ($fcg) = in-/^ {VuQ) (85) 

is simply related to the motion of the matter relative to the surfaces of constant 
temperature: 
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6.3. Adiabatic perturbations after decoupling 

After decoupling, in the free propagating domain, Rc(t) <^ 1, so the large-scale 
perturbation equations reduce to 



(rm) 



(86) 



S{rm) + H ( + 1 ) S(rm) = , 



(87) 



V{rm) + - 



V, 



1 1 



3 h ■('^"^ Aha^ 



(88) 



and in the matter dominated limit they can be integrated to give the following solutions 



A(^) = , (89) 

S(rm) =So + Bt'/' , (90) 
V^rm) = Ct . (91) 

These solutions demonstrate that after decoupling generic density perturbations do not 
satisfy the adiabatic condition S(^rm) = y{rm) = 0. This is due to the fact that the 
average velocity of the radiation does not proceed along geodesies, while the matter 
does. Thus any perturbation that starts off adiabatic at last scattering will not remain 
so. 

7. Conclusion 

In this paper we have calculated a fully covariant formula for the CMB temperature 
anisotropy improving on earlier work by Russ et al [^. This formulation has a number 
of distinct advantages over the more standard approaches as it is independent of gauge 
conditions, non-local splittings of spacetime, and related Fourier decompositions of 
perturbations around a FRW metric. Furthermore the results are relatively simple 
and easy to interpret. For scalar perturbations we recovered the dominant Sachs - 
Wolfe term, together with the Rees- Sciama effect which contributes to large scale CMB 
anisotropy only if the perturbations to the gravitational potential are non- stationary. 

We also examined the validity of the assumption that the density perturbations are 
adiabatic at decoupling and showed that if the surface of last scattering is correctly 
placed and the background is assumed to be a flat {K = 0) FRW model, then the 
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scalar (Sachs- Wolfe) contribution to large scale CMB anisotropics may be interpreted 
as being due to the motion of matter relative to the surfaces of constant temperature 
which define the surface of last scattering on scales where the instantaneous decoupling 
approximation applies. 
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Appendix A. Covariantly defined harmonics 



In the standard approach to cosmological perturbations [|l|, a harmonic decomposition 
of the perturbation variables is usually carried out using harmonics which are 
eigenfunctions of the Laplace - Beltrami operator on the three - hypersurfaces of constant 
curvature i.e. on the homogeneous spatial sections of FRW universes. In the covariant 
approach, the fluid four - velocity u"" is emphasized rather than an arbitrarily chosen 
spatial slicing, and quantities are defined by projecting orthogonal to -u" using the 
projection tensor hab- Covariant harmonics are therefore defined through operators 
constructed with the spatial (orthogonal to m") derivative ^^■'Va which are covariantly 
constant along the fluid flow lines (i.e. independent of proper time). In this section we 
will focus on scalar harmonics Q which are the eigenfunctions of the covariantly defined 
Laplace - Beltrami operator fl^ 



a? 



where A; is a comoving (i.e constant) eigenvalue. If z/ is a non- negative real wavenumber, 
then for a flat background {K = 0), it is associated with the physical wavelengths 
A = 27ra/z/, since in this case k = u, however for open models {K = —1) the spectrum 



of eigenvalues is given hj k = u + 1 |15 



The scalar harmonic Q can be used to define a vector 

Qa = -j'^'^yaQ (A2) 

k 

and a trace - free symmetric tensor 

Qab = ^^'V,(=^)V,Q + \KbQ . (A3) 

These harmonics are defined in order to have 

Q = Qa = Qab = 0, (A4) 
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so that they are covariantly constant along m". Finally, by taking the spatial divergence 
of ( |A3D the following relation between Qab and Qa is obtained 

a^^^V'Qab = -Ik-' {3K - k^) Qa . (A5) 

This is needed in section ^. 
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